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00 ; Abstract 

We study leptogenesis from the out-of-equilibrium decays of the lightest heavy 
C^l I neutrino N\ in the medium (low) temperature regime, T % 10 12 GeV (10 10 GeV), 

^ ' where the rates of processes mediated by the r (and /x) Yukawa coupling are non 

negligible, implying that the effects of lepton flavors must be taken into account. 
We find important quantitative and qualitative differences with respect to the 
case where flavor effects are ignored: (i) The cosmic baryon asymmetry can be 
enhanced by up to one order of magnitude; (ii) The sign of the asymmetry can be 
opposite to what one would predict from the sign of the total lepton asymmetry 



o 



^ ei; (iii) Successful leptogenesis is possible even with ei = 0. 

j§ 1 Introduction 



A very attractive mechanism for explaining the origin of the baryon asymmetry of 
the Universe (Yg = (ns — ng)/s — 8.7 x 1CT 11 ) is baryogenesis through leptogenesis 
P3 0j. Leptogenesis scenarios naturally appear within the standard model minimally 
extended to include the see-saw mechanism [3], because all the conditions [1] required 
for generating a cosmic lepton asymmetry are genericaliy satisfied in the decays of the 
see-saw related heavy singlet neutrinos: (i) The Majorana nature of their masses is a 
source of lepton number violation; (ii) Complex Yukawa couplings induce CP violation 
in the interference between the tree level and loop decay amplitudes; (iii) For a heavy 
Majorana mass scale ~ 10 11±3 GeV, sizable deviations from thermal equilibrium in the 
primeval expanding Universe can occur at the time the heavy neutrino decay. Partial 
conversion of the lepton asymmetry into a baryon asymmetry then proceeds by means 
of anomalous B + L-violating electroweak sphaleron interactions [5^ that are standard 
model processes. Qualitatively it is then almost unavoidable that a lepton (and hence 
a B — L) asymmetry is induced in the decays of the see-saw singlet neutrinos and, since 
no standard model reaction violates B — L, this asymmetry survives until the present 
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epoch. The question of whether leptogenesis is able to explain the puzzle of the baryon 
asymmetry of the Universe is then a quantitative one. 

The quantitative analysis of leptogenesis has become more and more sophisticated 
in recent years, taking into account many subtle but significant ingredients, such as 
various washout effects [HI 13 El 13 OH]- thermal corrections to particle masses and 
CP violating asymmetries [TTj . and spectator processes ^21 EI]- The latter are L 
and B — L conserving processes, such as standard model gauge interactions, some 
Yukawa interactions involving the heaviest fermions, and electroweak and strong non- 
perturbative 'sphaleron' interactions. These processes do not participate directly in 
the generation or washout of the asymmetries (hence the name 'spectator'), but have 
important effects in determining the asymmetric densities of the various particles, 
mainly by imposing certain relations among the chemical potentials of different particle 
species [12] ■ A detailed analysis of how spectator processes affect the washout back- 
reactions and concur to determine the final value of the baryon asymmetry has been 
recently presented in [T3] . 

In spite of all these refinements, one potentially very significant aspect of leptoge- 
nesis has been only rarely addressed fiJESEIlj and that is flavor. Neglecting flavor 
issues is, however, justified only if the process of leptogenesis is completed at a rather 
high temperature, T > 10 12 GeV. A reliable computation of the lepton and baryon 
asymmetries when leptogenesis occurs in the intermediate or low temperature windows 
must include flavor effects. In this paper, while taking into account all the effects 
discussed in jT2], we introduce in the analysis additional important phenomena that 
have to do with the flavor composition of the lepton states involved in leptogenesis. 1 In 
particular, we focus on the decoherence effects that are induced by the charged lepton 
Yukawa interactions on the lepton doublets produced in the decays of the heavy neu- 
trinos. As soon as these Yukawa interactions approach equilibrium, they act essentially 
as measuring devices that project all the lepton densities onto the flavor basis. Lepton 
number asymmetries and washout effects then become flavor dependent, and this can 
lead to a final baryon asymmetry that is different in size and possibly even sign from 
the one that would arise if flavor issues were irrelevant. 

The plan of this paper is as follows. In Section |2] we present the main physics 
ideas that underlie the most important effects of flavor in leptogenesis, and derive 
the relevant results in a qualitative way. In Section |H] we present the network of flavor 
dependent Boltzmann equations, including all the spectator processes discussed in [T3] . 
In Section 0] we separately analyze each of the relevant temperature regimes: we impose 
the appropriate equilibrium conditions and discuss their implications for the network 
of Boltzmann equations. We also present results for a set of representative flavor 
structures and for different temperature regimes. In Section |5] we compare our results 
to what is obtained when flavor effects are neglected or irrelevant, and we explain the 
main mechanism underlying the large enhancements of B — L that flavor effects can 
induce. In Section |H] we summarize our main results. 

our previous work |13| we imposed certain flavor alignment conditions, whereby the effects 
discussed in this paper become irrelevant. 
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2 The basic ideas 



In this section we study leptogenesis within temperature ranges well below 10 13 GeV, 
where lepton flavor issues can have an important impact on the way leptogenesis is 
realized. A proper treatment of the decay and scattering processes occurring in a 
thermal bath that consists of a statistical mixture of various flavor states should be 
carried out within a density matrix formalism, as discussed for example in Here 
we follow a simpler approach that is based on a physically intuitive formulation of the 
problem, and allows us to obtain all the qualitative features of the possible solutions. 
The results we derive here agree well with what will be obtained in Section 0]by solving 
numerically the detailed flavor- dependent Boltzmann equations that are derived in 
Section El 

2.1 Flavor CP violating effects 

The heavy singlet neutrinos N a needed for the see-saw model decay into lepton and 
antilepton doublets, implying lepton number violation. We denote the lepton doublets 
produced in the N a decays with the same index a = 1, 2, 3: 2 

r a = r(N a ^e a H), 

f Q ee T(N a ^F a H). 

CP- violation in N a decays can manifest itself in two different ways: 

1. Leptons and antileptons are produced at different rates, 

r a ^ f a . (i) 

2. The leptons and antileptons produced in N a decays are not CP conjugate states, 

CP(Q 5^ **• (2) 

If the rate of charged lepton Yukawa interactions is much slower than the rate of the 
heavy neutrino (N — £) Yukawa interactions, flavor issues can be neglected since, re- 
gardless of their flavor composition, £ a and £' a remain coherent states between two 
successive interactions. In this case only the effect in eq. (JTJ) is important. Indeed, 
leptogenesis studies have generally concentrated on this first effect. The situation is, 
however, different if leptogenesis occurs when the processes mediated by the r (and 
possibly fi) charged lepton Yukawa couplings are faster than the N — £ Yukawa interac- 
tions (a sufficient condition for this is that these processes occur at a rate comparable 
to the expansion rate of the Universe). Then, £ a and l' a are no longer the interacting 
states populating the thermal bath and £{ and ^ with i — r, (/i, e) should be consid- 
ered instead. The second effect in (J2J) can then become of major importance for the 
generation of cosmic asymmetries. 

2 The subindex in £1,2,3 bears no relation to the neutrino mass eigenstates ^1,2,3- 
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The amount of lepton asymmetry produced per N a decay (a = 1, 2, 3) is 

r — f 

- 1 - rv J- rv 



T a + T c 



(3) 



In order that non-vanishing e a asymmetry would arise, the two sets of three states t a 
and £' a cannot form an orthogonal basis because, to get CP violation from loops, it 
is required that the lepton doublet coupled to the external N a (the decaying heavy 
neutrino) couples also to a virtual N@ (with /3 ^ a) appearing in the loop. Assuming 
a hierarchical pattern for the heavy N a masses M a , and given the non-orthogonality 
of the i[ 2 3 states, the most natural situation is then that any asymmetry produced 
in A^2,3 decays is quickly erased by fast L violating interactions involving Ni. One can 
still envisage a situation in which £ 3 }fJL 2 is responsible for e 2 ^ 0, but an approximate 
orthogonality l 2 -L l\ prevents e 2 from being washed out by the Ni lepton number 
violating processes. In this case a lepton asymmetry produced in N 2 decays can survive 
until all the L violating processes involving Ni freeze-out. However, if this freeze- 
out occurs after lepton flavor dynamics has become important, then all the will 
be effectively projected onto £ T ,(u,e) flavor states. Then additional 'flavor alignment' 
conditions are necessary in order to preserve the e 2 asymmetry, as for example l\ oc £ e 
and ^2,3 -L ^e- Models that realize this scenario have been studied for example in 
[T4*l IT?! ITS] . In the following we disregard this possibility, and concentrate on the 
decay of the lightest heavy neutrino Ni (and on ei) as the dominant source of a cosmic 
lepton asymmetry. 

If during leptogenesis the (required) deviations from thermal equilibrium are not 
large, inverse decays and other washout process become important to determine the 
amount of e\ asymmetry that can be effectively converted into a baryon asymmetry by 
the B + L violating electroweak sphaleron processes. In the cases where the N\ heavy 
neutrinos have an initial thermal abundance, or when thermal abundance is reached due 
to inverse decays and other scattering processes, it is customary to express the present 
density of cosmic baryon asymmetry relative to the entropy density s as follows: 

— = -K s e 1 ri. (4) 

The numerical factor k s ~ 1.38 x 10~ 3 accounts for the B — L entropy dilution from 
the leptogenesis temperature down to the electroweak breaking scale, as well as for 
the electroweak sphalerons B — L — > B conversion factors. The factor 77, that can 
range between and 1, is the so called efficiency (or washout) factor and accounts for 
the fraction of lepton asymmetry surviving the washout processes. If N\ decays occur 
strongly out of equilibrium, then all back-reactions are negligible and 77 « 1. On the 
other hand, values of rj ~ 1CT 2 — 1CT 3 that are typical of strong washout regimes, when 
deviations from thermal equilibrium are mild, can still yield successful leptogenesis 
while ensuring at the same time that ub/s is largely independent of initial conditions. 

When Mi ^> 10 12 GeV, successful leptogenesis requires that the Ni — i\ couplings 
are sizable and, in particular, larger than all the charged lepton Yukawa couplings. 
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Then in the relevant temperature window, around T m Mi, charged lepton Yukawa 
processes are much slower than the processes involving N% (and also slower than the rate 
of the Universe expansion). In this regime, the composition of the two states £\ and £[ 
in terms of the lepton flavor states £i [i = r, fi, e) is irrelevant since the doublet states 
produced in the decays keep coherence between two different scatterings involving Ni. 
However, if leptogenesis occurs at lower temperatures, processes mediated by the r 
Yukawa coupling (and for T < 10 9 GeV also the processes mediated by the \i Yukawa 
coupling) become faster than the reactions involving N\. Then £\ and £[ lose their 
coherence between two subsequent L-violating interactions. Before they can rescatter 
in reactions involving N%, they are projected onto the lepton flavor basis with respective 
probabilities for each flavor i: 

K u = \(£i\£i)\ 2 and K li = \(? 1 \I i )\ 2 (CP(£ t ) = £ t ). (5) 

In the following, we drop the index 1 in Ku, Ku (as well as in the decay rates Ti, fi), 
leaving understood that we always refer to N± related quantities. We will also concen- 
trate exclusively on the medium and low temperature regimes, in which the charged 
lepton Yukawa interactions effectively 'measure', at least in part, the flavor composition 
of £±\ We distinguish the following possibilities: 

1. Alignment: if the i = r (or i — t, jj) Yukawa processes are in equilibrium, but 
both £\ and £' x are aligned or orthogonal to one flavor i (Ki — K^ — 1 or 0), 
leptogenesis reduces to a one-flavor problem, much like the unflavored case of the 
high temperature regimes. Different flavor alignments mainly affect just the way 
in which the lepton asymmetry gets distributed between lepton left- and right- 
handed states (the latter being sterile with respect to all L-violating processes) 
and this yields numerical effects not larger than a few tens of percent. A detailed 
analysis of the various aligned cases has been recently presented in |13j . 

2. Non-alignment, with only the r (but not the fi and e) Yukawa processes in equi- 
librium: the /i and e flavor components of £\ and £[ are not disentangled from 
each other during leptogenesis. It is then convenient to introduce two combi- 
nations of the /i and e flavors, 4 and £' b , that are orthogonal to, respectively, 
£i and £' x , thus implying Kf, = Kf, = 0. The two combinations of fi and e, £ a 
and £' a , that are orthogonal to, respectively, £b and £' b , satisfy K a = 1 — K T and 
K a = 1 — K T (in general CP(£' ab ) ^ £ a ,b)- In this case, even in the absence of 
any particular alignment condition, leptogenesis can be treated as an effectively 
two-flavor problem. 

3. Non-alignment, and both r and fi Yukawa interactions fast: all the lepton flavor 
components of £± are effectively resolved, and the thermal bath is populated 
by the (CP conjugate) flavor states £j and £-i (i = e,fi,r). The full three-flavor 
problem has to be considered in this case. (The onset of thermal equilibrium for 
the electron Yukawa coupling occurs at temperatures too low to be relevant for 
standard leptogenesis scenarios, but would not add qualitative changes to this 
picture). 
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In the rest of this paper we address the lepton flavor issues for leptogenesis in the 
above cases 2 and 3. We show how these issues have a significant impact on the way 
leptogenesis is realized in the intermediate and low temperature regimes. 

2.2 Lepton Flavor Asymmetries 

The CP asymmetry for N± decays into the £j lepton flavor is defined as 
■ = r(iV 1 -> £jH) - F(iV 1 -> IjH) Tj - Tj 

1 r+r r+r • I 

Since, by definition, Kj = Yj/Y and Kj = Yj/Y, we can conveniently express e\ as 

T Kj — Y Kj Kj + Kj Kj — Kj - AK 

A = ^ + f 3 = 3 2 3 ex + 3 2 1 ~ e x K] + (7) 

Here ei is the total asymmetry in iVi decays defined in eq. ©, AKj = Kj — Kj, and 
Kj = r°/r° represents the ratio of the (CP conserving) tree level decay amplitudes, 
with Kj = Kj. In the last equality in eq. (|ZJ), the first term, proportional to ei, corre- 
sponds to CP violating effects of the first type [eq. (JTJ)], the second term, proportional 
to AKj, is an effect of the second type [eq. (J2j)] since it vanishes when CP(^) = £±, 
and higher order terms of 0(AKj ■ ei) are neglected. 

For temperature regimes where the flavor states lj are the ones relevant for lepto- 
genesis, rather than the £ a states, eq. (@J should be replaced with 

n B Uf 

— = -k s J24vj, (8) 

S 3=1 

where rif = 2(3) is the relevant number of 'active' flavors in the medium (low) tem- 
perature regime. In these regimes, off-shell AL = 2 washout processes involving in 
particular ^2,3 are generally negligible, and hence all the relevant washout processes 
involve just the heavy neutrino N\ and are therefore associated with the flavor projec- 
tors Kj. Since, as will be discussed in more detail in Section the efficiency factors 
are, to a good approximation, inversely proportional to the washout rates, we can write 
r/j ~ min [r\jK Q -, l) , where r\ represents the washout factor one would obtain neglecting 
flavor effects. The value of r]j saturates to unity when Kj ~ r\ and, in the cases we 
have studied numerically, this typically occurs for Kj ~ few x 10~ 2 . In these cases, 
one can think of the decay of N\ — > £jH as one that proceeds much like strongly out- 
of-equilibrium decays. However, if Kj <C 1, the condition ^2 i K^ = 1 implies that for 
(one or two of the) other flavors, the back reactions are rather fast and in particular 
can be quite effective in populating the Ni states, so as to keep their abundance close 
to thermal during a relevant part of the leptogenesis era. As in the standard thermal 
leptogenesis scenarios, this ensures independence from initial conditions and, moreover, 
that a sizable amount of the asymmetry, much larger than a fraction Kj, can end up 
surviving in the £j flavor. Inserting eq. (JTJ) into eq. (JHJ) gives: 
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AKj 
i 2K, 



Kf > r] for all i 



(9) 



7]{n } - l)ei + K j e 1 + r)J2 i ^ j 



2Ki 



2 



We note the following: 

1. In the second line of this equation, that corresponds to the situation in which 
the value of rjj saturates to ~ 1, the first three terms are suppressed at least as 
1], while the last term is not. 

2. When K® <C rj, an alignment condition is approached [T!3] . Flavor effects become 
strongly suppressed here when all AKi vanish, as is always the case in the two 
flavor situations (see section HJ) . 

3. If AKi = for all i, that is in the absence of CP violating effects of the type of 
eq. (J2J), and if Kf > rj for all i, that is away from alignment conditions, then the 
final baryon asymmetry is always enhanced by a factor rif with respect to the 
case in which flavor effects are irrelevant (or are neglected). This enhancement 
occurs independently of the particular values of the K®, as is clearly seen from 
fig. HI and it holds also with respect to the aligned cases discussed in jTHj since 
alignment effectively enforces the condition n/ = 1. It is easy to understand the 
reason for the rif enhancement. When the total decay rates are projected onto the 
flavor i, each flavor asymmetry gets reduced compared to the total asymmetry: 
e\ ~ ei-Kj . However, this is compensated by a suppression of the corresponding 
washout factor, rji ~ f]/Kf. Then the sum over flavors yields the nf enhancement. 

4. When l x has approximately equal projections onto the different l± doublets, im- 
plying that Kf 1/nf for all i, the result of the previous point still holds, 
regardless of the particular values AKi 7^ 0. This can be easily seen from the 
first line of eq. by noting that Y2i ^K{ = 0. 

5. In the general situation, the terms proportional to AKi do not vanish. Let us 
stress that |Aifj/ei|, being the ratio of two (higher order) CP violating quantities, 
is not constrained to particularly small values and can well be sizable bigger than 
unity. This means that the effects of CP(£') ^ I can be the dominant ones in 
determining the size of the final baryon asymmetry. Moreover, since the signs of 
the AKi are not directly related to the sign of ei, it is clear that one cannot infer 
in a model independent way the sign of the cosmic baryon asymmetry only on 
the basis of e\. This situation (dominance of the AKj effect) is even more likely 
when one K® is very small and rjj ~ 1. Then, as can be seen from the second 
line in eq. Q, the term AKj/2 being not suppressed by the small value of rj can 
easily dominate over all other terms. 
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2.3 Dependence on Lagrangian parameters 

We now proceed to express the various relevant quantities - e\, e{, AKj and - in 
terms of the Lagrangian parameters. In the mass eigenbasis of the heavy neutrinos N a 
and of the charged leptons (e, denote the 577(2) singlet charged leptons), the leptonic 
Yukawa and mass terms read: 

C Y = ~M a N c a N c a - (A m N a U + h % e t t % & + h.c.) . (10) 

Explicit computation of the vertex and self-energy contributions to e\ with this La- 
grangian yields |TjJj: 

where xp = M^/Mf. The coefficient of (W^)pi within square brackets in eq. (jllj) comes 
from the leading term in the expansion for Xp ^> 1 of 



-^L + ^(l-(l + Xp ) In = -\x^ - \x^ + 0(x 

1 — xp \ xp J 2 p 6 



■5/2> 



(12) 



The first term in the l.h.s. originates from self-energy type of loop corrections, while the 
second term corresponds to the proper vertex correction. The coefficient of the second 
term (AA^)^ in square brackets in (fTTJ) comes from the self-energy type of diagram 
with 'inverted' direction of the fermion line in the loop, and corresponds to the leading 
term in the expansion (1 — xp)^ 1 = —x^ 1 + 0(x^ 2 ). Eq. (fTTj) holds when the mass 
splittings between the heavy neutrino masses are much larger than the decay widths, 
Mp — Mi 3> T^i. The lowest order expression for the flavor projectors reads: 

K _ ftO _ A lJ A lj / 13 x 
K 3 -K 3 - (AAf)ii . {16) 

By summing eq. (jllj) over the flavor index we obtain the total asymmetry: 

We do not need to separately calculate AKj in terms of the Lagrangian parameters 
since, from eq.((7|), AKj/2 — ej — ei Kj. As was first noted in ref. [JJJ], the term (XX^)ip 
contributes to the flavor asymmetry e\ in eq. (fTTj) . but - since (AA^) / gi(AA^)i / g is real 
- it does not contribute to the total asymmetry e\ in eq. (JHJ). We conclude therefore 
that it corresponds to a CP violating effect of the second type. 

In the basis that we are using, which is defined by diagonal matrices for the heavy 
neutrino masses M a and for the charged lepton Yukawa couplings hj, the X a j couplings 
can be conveniently written in the Casas- Ibarra parametrization [271] : 

1 



v 



MR^iU ] (15) 

J aj 
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where v = (H) is the Higgs vacuum expectation value, M = diag(Mi, M 2 , M 3 ) is the 
diagonal matrix of the heavy masses, m = dmg(m 1 ,m 2 ,m 3 ) is the diagonal matrix 
of the light neutrino masses, m x l 2 is an orthogonal complex matrix 

(R T ■ R = I) and U is the leptonic mixing matrix. With this parametrization, the term 
(AAt)^ that controls the total asymmetry in eq. (fH^ is given by 



(AA% = ^(E^%) > (16) 



while the term (A^-Ay) (AA^i, that gives the leading contribution to e{ in eq. (JTTJ), is 
given by 

(A/yAy (AAV = fer^iU ( E • ( 17 ) 



Eqs. and (J 17)) show that, if in the basis where the light and heavy neutrino 
Majorana masses are real, R were a real orthogonal matrix, two very intriguing conse- 
quences would follow: 

1. The total asymmetry would vanish, e\ = 0, while, in general, e{ ^ 0. This is a 
surprising possibility for leptogenesis. 

2. The asymmetries e{, which involve the imaginary part of eq. (|17|) . would depend 
only on the CP violating phases present in U, which (in principle) are measurable 
in low energy experiments. 

This situation would be quite different from the one usually considered in leptogenesis 
scenarios, in which the final baryon asymmetry, being proportional to e\, depends 
through the matrix R on all the low energy and high energy parameters. It remains to 
be seen if a real matrix R can naturally arise in some model. 



3 The network of Boltzmann equations 

We consider the scenario in which the heavy neutrino masses are hierarchical, Mi <C 
M 2) 3, and consequently the lepton asymmetry is generated mainly via the CP and 
lepton number violating decays of the lightest singlet neutrino N± to the lepton doublets 
l\ and £[. The important processes involving these states are the following (all the 7's 
below denote the thermally averaged rates): 

• N± decays and inverse decays, with rates 7^ = 7(^1 <-> £± H) and 7d = 7(^1 *-* 

Ah). 
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• AL = 1 Higgs- mediated scattering processes with rates such as 7s s = j(£i N\ <-> 
Q$t) and js t = 7(^1 Q3 *~* Nit), where Q 3 and t are respectively the third 
generation quark doublet and the top 577(2) singlet, as well as those involving 
gauge bosons, such as in £\Ni —>■ HA (with A = W 3,± or B). CP violating effects 
might be numerically important in these 2 2 scatterings jHJ E] . However, they 
do not add qualitatively new features to the analysis. Therefore we neglect CP 
violation in all these processes and accordingly we use the tree level expression 
for the flavor projectors Kj, Kj ~ Kj. 

• The s-channel scattering processes 7at s = j(£i H «-> £\ H) with on-shell N\ are 
already accounted for by decays and inverse decays. When the lepton doublets 
are projected onto the flavor basis, subtraction of these rates to avoid double 
counting must be carried out with care. 

• The off-shell scatterings 7^, involving the (pole subtracted) s-channel and the 
li-channel, as well as the t-channel scatterings jm = j(££ <-> H H), depend 
on all the projectors in a rather complicated way. However, since they are 
subdominant in the temperature ranges we are interested in, they can be safely 
neglected. 

In the temperature regimes in which the charged lepton Yukawa couplings become 
non-negligible (T <C 10 13 GeV), the corresponding interactions define a flavor basis 
for the lepton doublets. Then the decay rates and scattering processes involving the 
specific flavors £j and anti-flavors £i have to be considered, and the Boltzmann equations 
should track the evolution of all the relevant single-flavor asymmetries. To obtain the 
appropriate set of equations one can rely on the density operator approach discussed 
in ^3]. One defines a density matrix p as the difference between the density matrices 
for the leptons and for the antileptons, so that pa oc Y^., and normalizes p so that 
Yl/iPa — Yl- A very useful property of p is that the off-diagonal terms, with i ^ j 
(the coherences) vanish whenever the Yukawa interactions of one of the flavor states £{ 
or £j are in thermodynamic equilibrium, since this effectively projects the leptonic states 
present in the thermal bath onto those aligned or orthogonal to the flavors characterized 
by non negligible Yukawa interactions. 3 Another property that we use is that if the 
population of one state vanishes, pa = 0, then all the coherences pij associated with 
this state also vanish (this follows from e.g. the inequality paPjj > \Pij\ 2 )• These 
properties allow one to restrict the general equation for p to a subset of equations for 
the relevant flavor diagonal directions pa = Y^.. 

3 The transition region between the regimes where a specific lepton Yukawa coupling is completely 
negligible, and the one in which it mediates reactions in full thermal equilibrium, that is when Yukawa 
reactions for one specific flavor are approaching equilibrium, should be treated with care, since off 
diagonal entries in the density matrix might not be dumped fast enough to be neglected in the 
flavor dynamics |21j . However, within the temperature ranges in which the lepton Yukawa reactions 
for each flavor are fully in equilibrium or strongly out of equilibrium, our Boltzmann equations can 
be safely applied. 
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Following the approach outlined in [T5] . in writing down the Boltzmann equations 
we account for all the particle densities that are relevant to the washout processes. 
Moreover, in the evolution equations for the lepton flavor asymmetries we also include 
the term dYjf /dz that formally accounts for the fact that electroweak sphalerons 
constitute an additional source of lepton flavor violation. Then, for consistency, we 
also need to add the equation dYg/dz = dY^ w /dz to account for baryon number vi- 
olation by the sphaleron processes. Given that sphaleron interactions preserve the 
three charges Aj = B/3 — Li associated to anomaly- free currents, it follows that 
Yg /3 = Y^ w ■ By subtracting the equations for the lepton flavor densities from 
the equation for baryon number weighted by a suitable factor 1/3, we obtain the fol- 
lowing network of flavor dependent Boltzmann equations: 



dY, 



Ni 



dz 

dY Ai 
dz 



sHz 
-1 



Yn^ 

1 N 1 



1 ) ( 7D + 2 1Ss + H 



st) 



(18) 



Yi 



sHz 

2% + (yt - vq 



V e 1 



-1 I r 

Yn 1 



-(y £i +y H )K° 



Id 



+ 1 



Khst- 



Yi 



N- 



1 N l 



\vu +vt- vq 3 



K?lSs}, (19) 



where we have used the standard notation z = M\/T. In these equations, Y^ x = n^ 1 / s 
denotes the density of the lightest heavy neutrino (with two degrees of freedom) rel- 
ative to the entropy s, yx = (n x — Tlx)/ Tlx denote the asymmetries for the various 
relevant species X = £, H, t, Q% and all the asymmetries are normalized to the 
Maxwell-Boltzmann equilibrium densities. Notice that Y Ai is the Aj number density, 
also normalized to the entropy, with Yl. = (2y£ i + y ei )Y eq . The reaction rates are 
summed over initial and final state quantum numbers, including the gauge multiplici- 
ties. In the asymmetries yx, X = £, H or Q 3 label any of the two doublet components, 
not their sum, and hence we normalize yx to the equilibrium densities with just one 
degree of freedom. This is different from other conventions (e.g. the one used in |llj ) 
and allows us to keep the proportionality y x oc fix in terms of the chemical potentials, 
with the usual convention that e.g. /x^ is the chemical potential of each one of the 
two SU(2) components of the doublet £i. As already said, the subdominant AL = 2 
off-shell scatterings, and CP violation in (AL = 1) 2 <-> 2 processes, jHJ E] have been 
neglected in eq. (|19J1 . and we make two further simplifications: 

1. In eqs. (|18J ) - (|19JI and in what follows we ignore finite temperature corrections 
to the particle masses and couplings [TT] . In particular we take all equilibrium 
number densities equal to those of massless particles. 

2. We ignore scatterings involving gauge bosons, for whose rates no consensus has 
been achieved so far [TTJ |S] . They do not introduce qualitatively new effects and 
no further density asymmetries are associated to them. 



We would like to emphasize the following points regarding eq. fT9"j) : 



11 



1. The washout terms are controlled by the .fT-projectors. The sources of the asym- 
metry, on the other hand, receive additional contributions from the AK's. Note 
that AKi is not simply proportional to the respective Ki (see eq.flZJ)). This has 
important consequences: when, for at least one flavor, AKi /ex <L Kf, the results 
are qualitatively different from the case where, for all the flavors, AKi/ex <C Kf. 
In particular, the fact that the sign of AKi can be opposite to that of the Kf ex 
term opens up the possibility that leptonic asymmetry-densities of opposite sign 
are generated. 

2. Subtraction of the on-shell contributions from the s-channel Nx heavy neutrino 
exchange, that corresponds to AL = 2 s-channel scatterings, has to be performed 
with care. The cross sections are flavor dependent: r )^ i ~ 2 (iiH — > IiH) changes 
Li by two units, but other scatterings, 7^ i_1 (.^if — > Ylj&tjH) change Lj by 
only one unit. Furthermore, differently from the unflavored case, the AL = 
(ALi = 1) channels, ^ = \^H -> V . ^ijH), together with their asymmetries, 
must also be taken into account. Moreover, through inverse processes the asym- 
metries in the flavors j ^ i do affect the evolution of y^ (this is similar to the 
way uh, Vq 3 and y t act in the proper ALi = 1 channels in the second line of 
eq. (JUJ)). Nevertheless, we see from eq. (fT9|l that the result of the subtraction 
agrees with what one would obtain by naively generalizing from the case of the 
flavor independent equation (see e.g. ref. [TT]) to the flavor dependent case. 

3. Consider the case Ki = Ki = 1 (and hence K^ = Kj^i = 0). Then we have 
Ya-j* — 0. Since Yb-l = 2~^fc^A fe , the equation for Y"a i; coincides (as expected) 
with that for Yb-l in the unflavored case, or in the cases with flavor alignment, 
see eq. (13) in ref. [T3j . 

4. As a consequence of the sum-rules J2 i Ki = ^ i Ki = 1 (that imply £V AK { = 
and e\ = ex), the equation for Y B -l — J2k^k ^ as a simple form also in 
the general case Ki ^ Ki ^ 1: Similarly to the unflavored or aligned cases [T5] , 
this equation depends on ex as the source term, and on a single asymmetry- 
density jji = J2i^iV£i- The weighted sum of asymmetries, yt, represents the 
effective lepton doublet asymmetry coupled to the washout of Yb-l- Of course, 
all the complications related with the flavor structure are now hidden in yi whose 
detailed evolution is determined by the additional equations that still depend 
explicitly on Kf and e\. However, as we will see, there is always a point in the 
.fT-space for which yi oc Y B -l- In this particular situation, the equation for Y B -l 
decouples from the other equations and can provide a simple representative one- 
flavor approximation to the flavor dependent case, that still captures some of the 
main effects of flavor. 

The network of equations eq. (fT£j) and (fTTIj) can be solved after the densities y^ (or 
yi), yu and y t — yQ 3 are expressed in terms of the quantities Ya with the help of the 
equilibrium conditions imposed by the fast reactions, as described in the next section. 



12 



The value of B — L at the end of the leptogenesis era obtained by solving the Boltzmann 
equations remains subsequently unaffected until the present epoch. If electroweak 
sphalerons go out of equilibrium before the electroweak phase transition, the present 
baryon asymmetry is given, assuming a single Higgs doublet, by the relation [22] 

28 

n B = ^ u b~l- (20) 

If, instead, electroweak sphalerons remain in equilibrium until slightly after the elec- 
troweak phase transition (as would be the case if, as presently believed, the electroweak 
phase transition was not strongly first order) the final relation between B and B — L 
would be somewhat different |2Sj- 

4 The equilibrium conditions 

In this section we discuss the equilibrium conditions that hold in the different tem- 
perature regimes which can be relevant to study flavor effects in leptogenesis. Since 
leptogenesis takes place during the out of equilibrium decay of the lightest heavy right- 
handed neutrino Ni, i.e. at temperatures T ~ Mi, the relevant constraints that have 
to be imposed among the different particle densities depend essentially on the value 
of Mi. To allow for a straightforward estimate of the importance of flavor, we chose 
the relevant temperature windows as in ref. where flavor effects were irrelevant 
because of the imposition of alignment constraints, and also in presenting the results 
we follow closely that analysis. We use the equilibrium conditions specific of each 
temperature regime to express y^, yu and y t — yQ 3 in terms of the ^a/s. 

4.1 General considerations 

The number density asymmetries for the particles nx entering in eq. (fTT?|) are related 
to the corresponding chemical potentials through 

_ gxT 3 \ n x /T fermions, , . 

nx — n x — - | 2/jx/T bosons, (21) 

where gx is the number of degrees of freedom of X. For any given temperature regime 
the specific set of reactions that are in chemical equilibrium enforce algebraic relations 
between different chemical potentials [22J. In the entire range of temperatures relevant 
for leptogenesis, the interactions mediated by the top-quark Yukawa coupling h t , and 
by the 577(3) x SU{2) x U(l) gauge interactions, are always in equilibrium. Moreover, 
at the intermediate- low temperatures in which flavor effects can be important, strong 
QCD sphalerons are also in equilibrium. This situation has the following consequences: 

• Equilibration of the chemical potentials for the different quark colors is guaran- 
teed because the chemical potentials of the gluons vanish, \i g = 0. 
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• Equilibration of the chemical potentials for the two members of an SU(2) doublet 
is guaranteed by the vanishing, above the electroweak phase transition, of fi w + = 
—f^w- — 0. This condition was implicitly implemented in eq. (fT§j) where we 
used hq = jj UL = n dL , /i£ = /i eL = /i„ L and /j h = Hh+ = Hho to write the 
particle number asymmetries directly in terms of the number densities of the 
SU(2) doublets. 

• Hypercharge neutrality implies 

(^Qi + ~ ^* - M - ^0 + = , (22) 

i 

where Ui, d{ and e, denote the SU(2) singlet fermions of the z-th generation. 

• The equilibrium condition for the Yukawa interactions of the top-quark fi t = 
Ms + yields: 

yt-yQ 3 = Y' ( 23 ) 

where the factor 1/2 arises from the relative factor of 2 between the number 
asymmetry and chemical potential for the bosons, see eq. (J21|) . 

• Because of their larger rates, QCD sphalerons equilibration occurs at higher tem- 
peratures than for the corresponding electroweak processes, presumably around 
T s ~ 10 13 GeV |2*1*1 12*51 "2""]) and in any case long before equilibrium is reached 
for the r Yukawa processes. This implies the additional constraint 

^(2^-/^-/0 = 0. (24) 

i 

The condition in eq. (f2~""|) allows one to rewrite the r.h.s. of eq. ("T""]) in terms of only 
the two asymmetries ye i and yn- To express these asymmetries in terms of the Ya 4 we 
define two matrices C e and C H through the relations: 

y<, = -T, c ',^ »» = -£c»£t. (25) 

j j 

The matrices C and C H constitute a generalization to the case of flavor non-alignment 
of the coefficients q and c#, introduced in "U"]. The numerical values of their entries 
are determined by the constraints among the various chemical potentials enforced by 
the fast reactions that are in equilibrium in the temperature range (T ~ Mi) where the 
asymmetries are produced. The Boltzmann equations in eq. ()19|) can now be rewritten 
as follows: 



dF A , -1 (Yi 



N 



dz sHz 1 \ Y^ q 



5(c& + cf) 7b + 



l) (ct ~fss + 2p 7«J + (2 (% + Cf ) (j St + i 7 5.) 



yeq 



(26) 
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These equations are general enough to account for all the effects of the relevant spec- 
tator processes (Yukawa interactions, electroweak and QCD sphalerons) as well as for 
a general lepton flavor structure. 

The fact that the consequences of flavor cannot be easily read off from the system of 
coupled equations (j2Tlj) impedes a simple comparison with the results of the unflavored 
cases. We later show the results obtained by numerically integrating the set of coupled 
equations. But, in order to get some insight into the results, it is possible to introduce 
an approximation to the general equations (|26|) in the form of a one-flavor equation for 
Yb-l that accounts quite accurately for the numerical impact of flavor effects for two 
classes of models: 

1. Models in which N\ decays with approximately equal rates to all flavors (Kf ~ 
l/n,f for all 

2. Models in which all the flavor asymmetries e\ are dominated by the term Kf e\ 
[see eq. (J7J)]. 

From the discussion in section |2~21 it is clear that both kind of models have the common 
feature of being sensitive only to CP violating effects of the type in eq. (0). The 
approximation to be discussed below captures in full this type of effects, but it is blind 
to the CP violating effects of the second type in eq. ©, and therefore does not yield 
reliable results for the cases where the asymmetries are dominated by the effects of 
AK ^ 0. 

We proceed as follows. We consider particular flavor structures leading to 
values that satisfy the conditions ^ KfC^ = eg, independently of the value of j = 
1, . . .rif. We further introduce for the Higgs asymmetry an average coefficient ch = 
J2j Cf /rif and neglect the terms bCf = Cf — c#. (In all the cases that we consider 
SCj 1 /ch ~ 15%). Then, for this particular configuration, we can add the equations 
(J26|) over the flavors i to obtain an equation for Yb-l that is independent of flavor 
indexes: 



dY B - L -1 \(Y N \ Yb-l 

= 7w;\\Y*- l ) eilD+ T^ 



1 ,~ 

- (q + c H ) Id + 



1 J ( Q 75s + 1st) + (2 h + c H ) 1 75 t + ^7 5f 



yeg y~i ' 2 '°7 ' V " ' "° ' \ ' 2 



(27) 



Of course, it is useful to proceed in this way only to the extent that this special case 
is representative of a more general class of flavor configurations (corresponding e.g. to 
different sets of Kj). As we will see, this is indeed true for the two classes of models 
mentioned above. Since eq. (J27|) is similar in form to the one studied in ref. ^3]; where 
different coefficients q and c# corresponding to situations of flavor alignment were 
introduced, the values of eg and ch will give, by direct comparison, a measure of the 
possible impact of flavor effects in these models. 
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4.2 Specific temperature ranges and flavor structures 

We now discuss the flavor effects in the various relevant temperature ranges. In order to 
clearly show the impact of these effects, we conduct our analysis in a way that allows 
for a meaningful comparison with our results in [13J in which - due to appropriate 
alignment conditions - flavor issues played no role. Thus we use the same temperature 
windows (below T ~ 10 13 ) as in and (obviously) impose the same equilibrium 
conditions. 

As explained above, some of our main results can be understood from the specific, 
effectively one-flavor, cases that are described by eq. (J27j) . These cases are presented in 
Tableland can be easily compared with the corresponding results in Table 1 of [T3j . 
In the last column in Table Q we give the resulting B — L asymmetry. In the fourth 
and fifth columns, we give the values of eg and c#. We remind the reader that the sum 
Q + ch gives a crude scaling of the overall strength of the washout processes, while the 
ratio ch jet gives a rough estimate of the relative contribution of the Higgs asymmetry 
to the washout. 

For more general flavor configurations that do not belong to the two classes of 
models 1 and 2 above, and for which the approximation in eq. (J27j) does not hold, 
we present numerical results in a graphic form in Figs. ^ and EJ To disentangle the 
impact of the various effects from that of the input parameters, the B — L asymmetry 
is calculated in all cases with fixed values of rh\ = 0.06 eV and Mi = 10 11 GeV, where 
mi = v 2 (XX')n/ Mi determines the departure from equilibrium of the heavy neutrino 
Ni and provides an overall scale for the strength of the washout processes. The value 
rhi = 0.06 eV is intermediate between the regime in which departures from equilibrium 
are large and all washout effects are generally negligible (rhi < 10~ 3 eV) and the regime 
in which washout processes are so efficient that often the surviving baryon asymmetry is 
too small (rhi ~ 0.1 eV). The choice of this value is also motivated by the atmospheric 
neutrino mass-squared difference if neutrino masses are hierarchical. As concerns M 1; 
it is clear that the relevant temperature range is actually determined by it, yet we fix 
the value at Mi = 10 11 GeV in order to have a meaningful comparison of the various 
effects of interest. Namely, since in each regime considered the same asymmetries are 
produced in the decay of the heavy neutrinos, a comparison between the final values 
of B — L for the different cases can be directly interpreted in terms of suppressions 
or enhancements of the washout processes. We assume an initially vanishing value for 
Y/v and for all the particles density-asymmetries, but for rhi > 10 -2 eV the results are 
insensitive to the initial values. The values of the parameters adopted here as well as 
the initial conditions are the same as in [T3] . 

The four different temperature regimes that we consider are distinguished by the 
additional interactions that enter into equilibrium as the temperature of the thermal 
bath decreases. Of course, the most important of these reactions will be those mediated 
by the charged lepton Yukawa couplings. 

1) Bottom- and tau-Yukawa interactions in equilibrium (10 12 GeV ^ T < 10 13 GeV). 
Equilibrium for hb- and /^-interactions implies that the asymmetries in the SU(2) sin- 
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Equilibrium processes, constraints, coefficients and B — L asymmetry 
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Table 1: The temperature window is given in the first column, the relevant interactions 
in equilibrium in the second, and the constraints on the chemical potential in the third. 
(Chemical potentials are labeled with the same notation used for the fields: pQ i =Qi, 
/i^ =£i for the SU{2) doublets, p Ui = u i) /x^ = di, /i £i = ti for the singlets and = H 
for the Higgs.) We also give in the third column the K®- values that are compatible 
with eq. (J27|) . The values of the coefficients q and ch are given in, respectively, the 
fourth and fifth column while the resulting B — L asymmetry (in units of 10~ 5 x ei) 
obtained for rh\ = 0.06 eV and Mi = 10 11 GeV is given in the last column. 



glet b and e T degrees of freedom are populated. The corresponding chemical potentials 
obey the equilibrium constraints ^ = /ig 3 — and fi T = fi£ T — Possibly hb and h T 
Yukawa interactions enter into equilibrium at a similar temperature as the electroweak 
sphalerons [21] . However, since the rate of the non-perturbative processes is not well 
known, we first consider the possibility of a regime with only gauge, QCD sphaleron 
and the Yukawa interactions of the whole third family in equilibrium. Since electroweak 
sphalerons are not active, at this stage B = 0. 

The lepton (£ a ,£b,£ T ) and antilepton (£' a ,£' b ,£ T ) flavor bases are defined such that 
(£i\£ b ) = {£'i\£' b ) = 0, and hence \{£i\£ a )\ 2 = K a = l-K T and |(^K>| 2 = K a = l-K T . 
The 'prime' labeling the a and b antilepton doublets is a reminder of the fact that in 
general £' ab are not the CP conjugate states of £ a ^- Even though the charged lepton 
Yukawa interactions for both a and b states are negligible, the fact that Lb = pbv — 
ensures that the off-diagonal entries p a y and pb a > vanish. Given that Y^ b /Y eq = 2y^ b = 
0, the set of eqs. (ffilj) is reduced to just two relevant equations, and the equilibrium 
conditions (restricted to the a, r subspace) lead to the following values for Cf and C\y 

C H = 1(3, 4) and G l = - ( ^ ^ ^ . (28) 
16 V ' ; 32 V 1 12 / 
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The effective one-flavor approximation to non-alignment for this regime (the first 
line in Tabled)) corresponds to K® = 16/27, as follows from the condition £\ Ki Cfj = 
ce = a,r). This value ensures y oc Yb-l, giving in this case q = 2/9 (while 
ch = *Yuj =aT Cf /2 = 7/32). The enhancement of |Yb-l| in our representative case is 
of about 60% with respect to the aligned cases [T3] in the same temperature regime. 

In Fig. we display by a thick solid line the B — L asymmetry as a function of 
K®. This curve corresponds to K T = K T and is obtained by numerical integration 
of the system of equations (|26[). Three special cases are denoted with filled circles: 
The flavor-dependent case with = 16/27 (described in the Table) and the flavor- 
aligned cases K® = 0, 1 (discussed in [L3J). Two features that are expected from the 
qualitative discussion in section 12.21 are apparent in this figure: First, the effects of 
flavor misalignment are quite insensitive to the particular value of K® (as long as it is 
not too close to the alignment conditions). Second, the one-flavor approximation (the 
filled circle at K® = 16/27) provides a good estimate of Yb-l for generic values of K®. 

In the more general cases where K T ^ K T , several different parameters concur to 
determine the final value of B — L. One way to explore the possible results in this case 
would be to randomly sample the values of the Yukawa couplings and CP violating 
phases entering the expressions for K® a , AK ra and e x and present the results in the 
form of scatter plots. However, to gain a qualitative understanding of the various 
possibilities, we proceed in a simpler way. It is based on the fact that, if we take the 
tree-level Ni — > £ T H decay amplitude to zero (A ir — > 0) while keeping the total decay 
rates fixed [(AA^)n = const.], the r-flavor projector K® vanishes as the square of this 
amplitude (oc |Ai r | 2 ) while AK T vanishes as the amplitude (oc Ai T ). This suggests to 
adopt, as a convenient ansatz, the relation AK T oc yK^ also for finite values of K®, 
allowing us to explore the interplay between the contributions of AK T and of e\K® to 
e[ - see eq. (JJJ) - as well as to e\ = e\ — e[, by means of a simple two-dimensional plot. 
To this aim, we take AK T /2ei = k t ^J K® and we fix k t to the representative value of 
1/4. Thus, for the practical purpose of carrying out a qualitative survey of the possible 
different situations, we adopt the following ansatz: 

ansatz : ^- = K° t + -Jk*. (29) 

With respect to the hierarchy between the two different CP violating effects of eqs. (JTJ) 
and (J2J), this simple relation has the nice property of covering all the interesting possi- 
bilities: 

1. For K® < 1/16, the two terms K® e\ and AK T are comparable in size. 

2. For 1/16 ^ K® < 3/4, K®e\ and K®e\ dominate the respective asymmetries. 

3. For K® > 3/4, we enter a regime in which AK T /2 ^> K®ei and, since AK a = 
—AK T , the asymmetry e", that is largely dominated by the AK a term, has now 
the opposite sign with respect to e x and e[. 
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Figure 1: The value of the final I^b-lI (i n units of 10~ 5 |ei|) as a function of K® in 
the regimes 1 (solid lines) and 2 (dashed lines), computed with Mi = 10 11 GeV and 
rhi — 0.06 eV. The filled circles (regime 1) and squares (regime 2) give the aligned 
cases (K® = 0, 1) discussed in ref. [T3] and the non-aligned cases given in the first two 
rows in Tabled The thick lines correspond to non-aligned models for which AK T = 0, 
implying t\/ti = K®. The thin lines give an example of the results for K T ^ K T 
assuming AK T /2ti = ^JK^/A. The corresponding values of e[/ei are marked on the 
upper x-axis. The arrows with labels (a),(6),(c) and (d) correspond to the four panels 
in figure (j2J). Note that Yb-l changes sign in (c). 



In Fig. ^ we show with the thin continuous line the results for this more general 
case. On the upper x-axis we have marked for reference the value of e[/e% corresponding 
to each different value of K®. The most peculiar features are the two narrow regions 
marked with (b) and (d) where is strongly enhanced. The enhancement takes 

place at values of not far from alignment. In particular, the steep rise of \Yb-l/^x\ 
close to K® ?a 1 can reach values up to one order of magnitude larger than the vertical 
scale of the figure. (Note, however, that the ansatz eq. (}2l3j) does not yield the required 
behavior AK a = —AK T = for = 1 — K® = 0, and therefore the thin lines 
in the plot should not be extrapolated to K T ~ 1.) It is worth noticing that that 
the two peaks correspond to values of Yb-l I ^-i of opposite sign, since the asymmetry 
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changes sign in (c). Qualitatively similar effects occur in the lower temperature regimes 
discussed below, and therefore we postpone the analysis of these results to section 
As concerns the quality of our model independent approximation, we see from Fig. ^ 
that, in this more general case, it can provide a reasonable estimate of Yb-l only for 

~ K® ~ 1/2, that is when Aq decays with similar rates into l T and l a . 

Given that the final value of Yb-l depends linearly on e[ (see eq. ISJ), we see from 
eq. (J7J) that, for a fixed value of K®, Yb-l is linear in AK T . Therefore, from the 
two thick and thin lines in Fig. ^ that respectively give Yb-l / ^i for AK T = and 
for AK T /(2ei) = k t a/A"°, one can easily infer, for each A"°, the value of Yb-l/ci 
corresponding to any other value of AK T . In particular, by reflecting the thin line with 
respect to the thick one, one can figure out the results one would obtain for k t < 0. 

2) Electroweak sphalerons in equilibrium (10 11 GeV ^ T < 10 12 GeV). 

The electroweak sphaleron processes take place at a rate per unit volume Y /V oc 
T 4 atyy \og(l/a>w) j27[ I28| I29j. and are expected to be in equilibrium from tempera- 
tures of about ~ 10 12 GeV, down to the electroweak scale or below [21]. Electroweak 
sphalerons equilibration implies 

5^(3/i Qi + /i A ) = 0. (30) 

i 

As concerns lepton number, each electroweak sphaleron transition creates all the dou- 
blets of the three generations, implying that individual lepton flavor numbers are no 
longer conserved, regardless of the particular direction in flavor space along which the 
doublet l\ and £[ lie. While the electroweak sphalerons induce L b ^ 0, the condition 
A b = is not violated, and hence eqs. (pjj) still consist of just two equations for Y& a and 
for Y& T . As concerns baryon number, electroweak sphalerons are the only source of B 
violation, implying that baryon number is equally distributed among the three quark 
generations, that is B/3 in each generation. This modifies the detailed equilibrium 
conditions for the quark chemical potentials. Besides these differences, the analysis 
follows closely the one carried out in the previous regime. The coefficients Cf and CL 
are given by 

C H = J_(41, 56) and C l = — ( ^ . (31) 

230 V ; 460 V - 9 156 / 

In Tabled w e give the values of q and c# that correspond to the approximation in 
eq. (J27|) . For the models described by this approximation, flavor misalignment can 
induce an O(80%) enhancement of the final B — L asymmetry compared to the aligned 
cases discussed in jT^]. The source of this enhancement is mainly the suppression by 
a factor ~ n/ = 2 of the washout processes that is apparent in the reduced value of 
C£. In Fig. ^we give the final value of as a function of K® for AK T = (thick 

dashed line) and for AK T ^ (thin dashed line). The qualitative dependence of Yb-l 
on is quite similar to regime 1: If the only CP violating effects are of the type 
in eq. ((TJ) (CP(^) = £i), the results are approximately independent of the particular 
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value of K®, while relaxing this condition results again in a strong dependence on K® 
and possibly strong enhancements. This general case is again well approximated by 
the values of q and ch given in the Table when K® ~ K® ~ 1/2. 

3) Second generation Yukawa interactions in equilibrium (10 8 GeV ^ T < 10 11 GeV). 

In this regime, h c , h s and, most importantly, Yukawa interactions enter into equi- 
librium. Given that the electron remains the only lepton with a negligible Yukawa 
coupling, the Yukawa interactions completely define the flavor basis for the leptons as 
well as for the antileptons (that are now the CP conjugate states of the leptons). Cor- 
respondingly, the lepton asymmetries are also completely defined in the flavor basis. 
In this regime, the coefficients and Cf projecting the asymmetries (yi e , yi , yi T ) 
and y H onto (Y Ae , Y A/i , Y At , ) are: 

j 1 / 906 -120 

C = -(37, 52, 52) and C = — ^ -75 688 -28 | . (32) 

For the approximation of eq. (J27|) . the coefficient c~i in Tabled is reduced by at least a 
factor of three with respect to the values of q in the aligned cases analyzed in [T3] , and 
the final value of Yb-l gets correspondingly enhanced. This is precisely the rif ~ 3 
enhancement expected from the qualitative discussion in section l2~2l The estimate of 
Yb-l based on the one-flavor approximation is, again, rather precise for the class of 
models for which AKi = 0, independently of the particular values of K® T (as long 
as they are not too close to or 1, where approximate alignment could induce effects 
that suppress Yb-l)- I n the more general case with AK e> ^ T ^ 0, the approximation 
is again reliable in the region of approximately equal flavor composition for l\ and £[, 
that is around K® ~ K® ~ ~ 1/3. As regards regions in K space away from equal 
flavor composition, we expect that the enhancements observed in the rif = 2 cases, 
that are especially large when one Kf is small (see fig. HJ), could be even larger in this 
case since now two flavor projectors can be simultaneously small. 

4) All SM Yukawa interactions and electroweak sphalerons in equilibrium (T 10 8 GeV). 

Equilibration of the Yukawa processes for all quarks and leptons, including the electron, 
occurs only at temperatures T < 10 6 GeV, which are too low to be relevant for leptoge- 
nesis (at least in the standard scenarios). We nevertheless analyze this regime, mainly 
for the purpose of comparison with analysis that assume that all Yukawa interactions 
are in equilibrium, but do not include the flavor effects [T2J HHj . 
The coefficients C H and C e are given by 

C H = ^(1, 1, 1) and C l = ^\ -16 221 -16 | . (33) 

Given the symmetric situation of having all Yukawa interactions in equilibrium, the 
point of equal flavor composition for the £p doublets (K e = = K T = 1/3) defines 
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a fully flavor-symmetric situation for which Y& T = Y& = Y& e = Yb-l/ 3. It is then 
straightforward to see that this point corresponds to the condition y e = £\ Kfy^ = 
—C£YB-i/y eq that defines our one-flavor approximation. The symmetric situation im- 
plies that an exact proportionality is also obtained for yu = —chYb-l/Y^ (indepen- 
dent of the particular values of K® ). In agreement with the qualitative analysis in 
section |2~21 flavor effects suppress q by the large factor 3.5 (~ rif) compared to the 
corresponding aligned case ^3]. in the cases restricted by AKi = for which the 
results of the one-flavor approximation hold, the final value of Yb-l is enhanced by a 
factor ~ 2.5 with respect to the aligned case jTSj. Again, much larger enhancements 
are possible if the condition Ki = Ki is relaxed. 

5 Discussion 

In this section we explain some generic features of our results. We refer here to results 
that, while depending on the specific flavor structures, are qualitatively similar in all 
the temperature regimes in which flavor effects are important. For the purposes of 
this discussion, it is important the fact that in all the regimes considered the final 
asymmetries are inversely proportional to the rates of the corresponding washout 
processes. This can be demonstrated along lines similar to those given in Appendix 2 
of ref. Jl] for Yb-l- Note that we have already used this result extensively for the 
qualitative discussion in section l2~2l In all the temperature regimes where flavor effects 
are important (Mi <C 10 13 GeV), the washout rate having the strongest impact on 
the final value of the asymmetries is 7^, that is the on-shell piece of the AL = 2 
scattering (the term proportional to Kf'jB in the first line of eq. ([19))) that has a 
Boltzmann suppression factor exp(— z), similar to the AL = 1 rates. Hence, the proof 
given in ref. [Hj for the case of AL = 1 washout dominance and small departure from 
equilibrium, holds also for the cases considered here, and applies to the single flavor 
density asymmetries as well. 

As is written explicitly in eq. (JBJ), the final values of the Y^ are determined by 
the asymmetry parameters e\ and by the washout factors rji. The washout factors 
are related to the various lepton number violating processes of eq. aud hence 

depend on the factors Kf that control the overall strength of the washouts, as well 
as on the matrix coefficients C e and C H defined in (125)) . However, the matrices C 
in eqs. PS]). (f3*T|). (f3l2|) and (jBlIj) are approximately diagonal and the diagonal terms 
in each matrix are not very different from each other (actually, for each temperature 
regime, Cf { equals the value of q in the case that Kf = 1, see |H|). If we make 
the crude approximations of (i) neglecting the off-diagonal elements, and (ii) taking 
the diagonal elements equal (that is C e oc 1), and we also note that the same Higgs 
asymmetry yn enters the equation for the different flavors, we are led to conclude that 
the relative values of the rji are determined mainly by the values of the Kf. More 
precisely, given that the amount of flavor-asymmetries surviving the washout are, to a 
good approximation, inversely proportional to the washout rates, the washout factors 
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Figure 2: Y^ T je\ (short-dashed line), Y/\ a jt\ (dotted line) and Yb-lAi (solid line) as 
a function of z = M\/T computed with Mi = 10 11 GeV and rh\ = 0.06 eV, in the 
temperature regime 2. The panels refer to the labels in fig. 1. The examples include 
both a constrained case, AK T = (a), and general cases with AK T /2e = a/ 'K^/A 
(b,c,d). The K° T value is (a) 0.1, (b) 0.02, (c) 0.84 and (d) 0.98. Note that the vertical 
scale in (d) is doubled with respect to the other three cases. 

obey the approximate proportionality rji oc l/K{. This results constitutes the basis of 
the approximate expressions for ris/s in eq. that allowed us to estimate qualitatively 
all the most important effects of flavor. We rely again on this approximation in the 
following discussion of the effects of different flavor structures and, in particular, in 
our attempt to understand the different behaviors that are apparent from the points 
labeled (a), (6), (c) and (d) in Fig. □ 

(a) This case is defined by the condition Ki = or, more generally, by the con- 
dition AKi €i K® for all relevant i. As we have seen, in this case the final Yb-l 
shows a sizable enhancement with respect to the aligned cases in the same temperature 
regime, and its value is approximately independent of Kf. This behavior corresponds 
to the ~ rif enhancement discussed in section l2~2l since Ki = Ki implies e\ = Kf e, 
it follows that, to first approximation, each Y&. oc e\r]i (i = a,r in regimes (1) and 
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(2), i = e,fi,r in regimes (3) and (4)) is independent of Kf and, moreover, that the 
final values of the YAr as y mme ^ r i es are * n general quite similar. The enhancement of 
Yb-l by a factor ~ rif with respect to the aligned cases [T3] is indeed confirmed by the 
numerical results given in Table E](w/ — 2 in regimes (1) and (2), rif = 3 in regimes (3) 
and (4)). This case is illustrated in more detail in FigJ2]a. The figure shows the evolu- 
tion of Y"A a and Y& T with z in the regime 2, for K® = 0.1. It is apparent that in spite of 
the large difference in the values of K® and K®, that result in quite different evolutions 
for the two asymmetries, their final values are very similar and approximately equal to 
Yb-l/2. 

Note that this result relies in a crucial way on the assumption that the source 
terms e\ and the washout rates are both proportional to K? . If this condition is not 
realized, as in the general cases when K ^ K, more subtle effects become important 
in determining the final values of the asymmetries Y"a; and, as is apparent from Fig^ 
the final value of Yb-l acquires a rather strong dependence on the values of Kf. Let 
us now discuss a few examples of this more general case. 

(b) In this case, AK T /2 and e\ K2 are of the same order and of the same sign. 
The details of this case are illustrated in FigJ2]& where we present the behavior of Ya t 
and Y Aa for K° = 0.02, that corresponds to AK T /2e = ^K^/A ~ 0.035. The picture 
explains rather well the origin of the first peak, labeled with (b) in Fig^ We see that a 
deviation at the level of just a few percent from exact alignment (K° = 0) is sufficient 
to start populating the asymmetry in the r flavor even if the rate, which is suppressed 
by K® + a/ K®/4 ~ 1/20, is rather small. For the first part of the leptogenesis era, 
Ya t gives only a minor contribution to the total B — L asymmetry. However, while 
the washout in the direction £ a (-L £ T ), that is controlled by K® « 1, proceeds with full 
strength, there is practically no washout for the r component and, eventually, while 
Ya b ends up being almost completely washed out, it is Y"a t that determines the final 
value of B — L at the end of the leptogenesis era. 

(c) This intriguing case corresponds to e[r] T ~ — e\r\ a : The final values of the two 
leptonic asymmetries are approximately equal in magnitude but of opposite signs. 
This implies that Yb-l can vanish even when lepton flavor asymmetries are sizable. 
Eq. (JHJ) helps us to understand this case in a simple way. Adopting the approximation 
r/i ~ 7)/K$ discussed above we have: 



2 6!- 



1 \ AK T 



1-K2 K2 



(34) 



With our particular ansatz, AK T /2ei = ^K2/4:, the r.h.s. of eq. ()34j) vanishes for 
K2 ~ 0.9, in qualitative agreement with the point labeled with (c) in Fig. ^ The 
evolution with z of the flavor asymmetries is illustrated in FigJ2]c. During most of the 
leptogenesis era, Yb-l is rather large. However, in the end its value is driven to zero, 
even if the final values of the two lepton asymmetries remain quite sizable (about twice 
larger than case (a)). For larger values of K® the total asymmetry Yb-l crosses zero 
and changes sign, and we enter a different regime to be discussed next. 
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(d) In this example we have AK a /2ei ~ -0.25 and = 0.02, so that AK a /2e 1 
is much larger in absolute value than K®, and is negative. This situation means that 
e" has a sign opposite to both e\ and €\. This case is illustrated in Fig. in which 
the vertical scale is doubled with respect to the previous three cases. Initially, a large 
density-asymmetry starts being built in the r flavor, due to the large value el/ei ~ 1.23. 
However, its growth is kept under control by the washout effects that, with K® « 1, 
are unsuppressed. In contrast, for Y^ a the washout effects are strongly suppressed 
by K® ~ 0.02 and eventually this asymmetry largely prevails, driving Yb-l to values 
several times larger than in all the previous cases, and having opposite sign. Notice 
that in this case the difference in sign between Y^ a and Ya t gives large cancellations 
between the respective contributions to the asymmetry y^, and this implies that the 
Higgs washout effects are accordingly suppressed. 

6 Conclusions 

We have shown that the effects of flavor can have dramatic consequences in leptogenesis 
scenarios. This occurs due to the way in which charged lepton Yukawa interactions in 
thermal equilibrium affect the flavor composition of the leptonic density-asymmetries, 
that determine the washout processes. For these effects to be significant, at least one 
leptonic Yukawa interaction (h T ) must be in equilibrium. This happens if Mi, the 
lightest heavy neutrino mass that determines the temperature at which leptogenesis 
takes place, is light enough: Mi < 10 13 GeV. Moreover, for the flavor effects to have 
an impact, one needs to be in the strong washout regime: rh\ > 2 x 10 -3 eV. 

The main consequence of the flavor effects is that, in generic flavor non-aligned 
cases with strong washouts, the final asymmetry is typically enhanced by a factor nf. 
We have rif = 2 when just h T is in equilibrium, i.e. for 10 9 GeV < Mi < 10 13 GeV, 
and nf = 3 when also is in equilibrium, i.e. for Mi < 10 9 GeV. 

In addition to the total asymmetry €\ associated to N\ decays, the individual asym- 
metries e{ can play a crucial role. These flavor asymmetries introduce a qualitatively 
new ingredient, AAj, which is the contribution associated to the fact that the leptons 
and antileptons produced in Aq decays are, in general, not CP conjugates of each 
other. When these new contributions are non- negligible, much larger enhancements 
become possible, especially when at least one lepton flavor £j is weakly coupled to the 
decaying Aq (typically this means K® ~ 77 <§C 1) . Since the signs of AKj and of e± need 
not be the same, the final asymmetry can have a sign unrelated to that of e±. Actually, 
successful leptogenesis is possible even with e\ = 0. Scenarios in which e x = while 
e{ 7^ entail the possibility that the phases in the light neutrino mixing matrix U are 
the only source of CP violation. 
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